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ABSTRACT 

We report on the discovery of an infrared cometary nebula around PSR J1549— 4848 in our Spitzer 
survey of a few middle-aged radio pulsars. Following the discovery, multi-wavelength imaging and 
spectroscopic observations of the nebula were carried out. We detected the nebula in Spitzer IRAC 
8.0, MIPS 24 and 70 /Ltm imaging and in Spitzer IRS 7.5-14.4 fxm spectroscopic observations, and 
also in the WISE all-sky survey at 12 and 22 /zm. These data were analyzed in detail, and we find 
that the nebula can be described with a standard bow-shock shape, and that its spectrum contains 
polycyclic aromatic hydrocarbon and H2 emission features. However, it is not certain which object 
drives the nebula. We analyze the field stars and conclude that none of them can be the associated 
object because stars with a strong wind or mass ejection that usually produce bow shocks are much 
brighter than the field stars. The pulsar is approximately 15" away from the region in which the 
associated object is expected to be located. In order to resolve the discrepancy, we suggest that a 
highly coUimatcd wind could be emitted from the pulsar and produce the bow shock. X-ray imaging 
to detect the interaction of the wind with the ambient medium and high-spatial resolution radio 
imaging to determine the proper motion of the pulsar should be carried out, which will help verify 
the association of the pulsar with the bow shock nebula. 

Subject headings: infrared: ISM — ISM: structure — stars: individual (PSR J1549— 4848) — stars: 
neutron 



1. INTRODUCTION 

Interstellar shocks are seen in association with var- 
ious a strophysical objects, f rom wind-blowing massive 
stars (Ivan Buren et al. Ml990[ ). mass-ejecting giant stars 
(e.g., [Martin ct al. 200 7[), expandi ng supernova remnants 

\ (SNRs; e.g., iM orton e t al.l I2007D . to energetic pulsars 
(jGaensler fc Slana j2006[ ) . In these objects, fast moving 
ejecta collide with the surrounding interstellar medium 

' (ISM), and the strong interaction drives a shock front. 
Depending on the kinetic energy carried into the shock 
and local ISM conditions, the postshock hot material 
can be bright at different wavelengths. Observations of 
shocks not only help study shock physics, but also pro- 
vide probes of the associated objects and ISM. 

Depending on temperature T, the sound speed Cs in 
the ISM [cs ~ 10(r/10'' K)i/2 km s"^] is in the range of 
1-100 km s^"'^. As a result, when stellar objects have high 
space velocities V-t , they can be moving supersonically in 
the ISM. In such cases, a shock will appear as a cometary 
structure (so-called bow shock), as interste llar material is 
swept up into a dense cone- type shell (e.g., 'Wilkin 1996). 

One type of bow shock is known to be associated with 
pulsars. Haying s pace velocities of >100 km s~^ (e.g., 
iHobbs et al.ll2005h and strong winds (which carry much 

of the spin-down energy E of pulsars), pulsars moving 
in the ISM are candidates to drive bow shocks. Thus 

^ Shanghai Astronomical Observatory, Chinese Academy of 
Sciences, 80 Nandan Road, Shanghai 200030, China 

^ Physics Department, University of WisconsinMilwaukee, 
Milwaukee, WI 53211, USA 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA 02138, USA 

* Las Campanas Observatory, Observatories of the Carnegie 
Institution of Washington, La Serena, Chile 

^ Department of Physics, McGill University, 3600 University 
Street, Montreal, QC H3A 2T8, Canada 



far, six pulsars have been found to have a typical bow 
shoc k, revealed by the shock s' optical Balmer line emis- 
sion (jGaensler fc Slanjl2006[ ). The detected line nebulae 
are understood to arise from de-excitations of neutral 
H atoms in the shocked ambient gas, following the pro- 
cesses of collisional excitation or charge-exchange (e.g., 
Bucciantini & Bandiera 2001). In addition, a pulsar bow 
shock may also be revealed by the detection of a termina- 
tion shock of the pulsar wind. The termination shock lies 
inside a bow shock, and can appear bright at X-ray or 
radio ener gies due to synchro tron radiation (|Kaspi et al.l 
2006; Gae nsler fc Sla^ [2006h . 

In our Spitzer Infrared Array Camera (IRAC) survey 
of 7 relatively young pulsars, for the purpose of prob- 
ing the general existence of debris disks around isolated 
neutron stars (Wang et al. 2013, in preparation), a mid- 
infrared (MIR) cometary nebula was serendipitously de- 
tected near one of the pulsars, PSR J1549— 4848 (here- 
after J1549). Naively thinking that the MIR structure is 
associated with the pulsar, we would have found the first 
MIR bow shock that arises from the interaction between 
the pulsar wind and ISM. We thus made follow-up 24 
and 70 /xm imaging and 7.5-14.4 /im spectroscopic ob- 
servations with Spitzer and optical observations with the 
Magellan telescopes after the initial detection of the neb- 
ula at 8.0 /im. However, our detailed data analyses raise 
questions about the association between the cometary 
nebula and the pulsar. On the basis of our observational 
results, the origin of the nebula is not certain. Here we 
report our multi- wavelength observations, data analyses, 
and results. The observations and data reduction are 
described in § [2l and the data analyses and results are 
presented in §[3]and §31 respectively. We discuss the pos- 
sible origin of the cometary nebula in §[5] and summarize 
the results in § [6l 
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Fig. 1. — Spitzer IRAC 8.0 fim (left panel) and MIPS 24 (third panel) and 70 fira (right panel) imaging of the cometary nebula. The 
WISE third-band 12 /im detection of the nebula is shown in the second panel. In the left panel, the position of J1549 is marked by a cross 
sign, and one nod position of the slit (which has a size of 57"x3'.'7) in our Spitzer IRS observations is indicated with a box region. 



2. OBSERVATIONS AND DATA REDUCTION 

Our Spitzer and ground-based observations of the 
cometary nebula are summarized in Table [1] Below we 
describe them and the related data reductions. 

2.1. Spitzer Imaging 
2.1.1. IRAC And MIPS Imaging 

As part of our survey program, the initial Spitzer ob- 
servations of J 1549 were carried out on 2007 September 
12. The imaging in strument used wa s the Infrared Ar- 
ray Camera (IRAC; iFazio et all 120041 ). The field of our 
target was imaged in the pair of 4.5 and 8.0 /im chan- 
nels of the IRAC simultaneously. The detectors at the 
short and long wavelengths were InSb and Si:As devices, 
respectively, with 256x256 pixels and a plate scale of 
l'.'2/pixel. The field of view (FOV) was 5'.2x5'.2. The 
frame time was 30 s, with 26.8 effective exposure time 
per frame. The total exposure time at each channel was 
17.87 min. A cometary nebula was detected in the 8.0 
/xm image (Figure [T]) around J1549, but not in the 4.5 
fim image. 

To study the nebula and derive its properties, we re- 
quested Spitzer observations at the longer wavelengths. 
Using Multiband Imaging Photometer for Spitzer (MIPS; 
iRieke et al.|[200l . 24 and 70 /im broad-band imaging of 
the source field was carried out on 2008 August 30. The 
24 /xm detector was a 128x128 pixel Si: As array, with a 
pixel size of 2'.'55 and a 5.'4x5f4 FOV. The 70 /im detec- 
tor was a 32x32 pixel Ge:Ga array. We chose its wide 
FOV observing mode for imaging, which provided a pixel 
scale of 9'.'85xl0'.'06 and a 5'.2x2'.6 FOV. At 24 and 70 
/xm, the exposure times were 11.8 and 10.5 min, respec- 
tively (frame times ~ 3 sec). 

2.1.2. Data Reduction 

We started with the Basic Calibrated Data (BCD), 
provided by Spitzer Science Center (SSC). The data were 
produced from raw images through the IRAC and MIPS 
data pipelines (IRAC software version was S16.1.10 and 
MIPS software version was S18.1.0) in SSC. The detailed 
reduction in the pipelines can be found in the IRAC Data 



Handbook (version 3.0) and MIPS data Handbook (ver- 
sion 3.3.1). 

For IRAC data, we first used the IRAC artifact miti- 
gation program, distributed as Contributed Software by 
SSC, to clean the BCD images. Several bright stars in 
the target field caused column pull-down, row pull-up, or 
optical banding artifacts. We then corrected each BCD 
image for array location dependence, using the array lo- 
cation correction images provided by SSC. To combine 
BCD images into a final post-BCD (PBCD) mosaic at 
each channel, we used MOPEX, an SSC's package for re- 
ducing and analyzing imaging data. 

For MIPS 24 /im data, the BCD images do not have 
any artifacts that need corrections. We used MOPEX to 
combine the BCD images into a final PBCD mosaic. The 
MIPS 70 /(m BCD images contain horizontal and ver- 
tical strips, which were due to response variations and 
stim latents (stims were MIPS's built-in calibration light 
sources, which flashed during an observation sequence; 
for details, see MIPS Data Handbook), respectively. We 
used the GeRT package, which is provided for MIPS 70 
& 160 /xm data reductions, to clean up the BCD images. 
During the reduction, the region that contains our source 
was masked out. The BCD images were combined into 
one final PBCD mosaic using MOPEX. 

2.2. Spitzer Spectroscopy 
2.2.1. IRS Observations 

Our Spitzer spectroscopic observations were made on 
2009 September 14. The instrument was the Infrared 
Spectrograph (IRS; IHouck et al] I2004D . We used the 
Short-Low (SL) module of the IRS for the observations, 
which provided a wavelength coverage of 7.5-14.4 /im 
(first order) and a resolving power of 61-120. The de- 
tector was a 128x128 arsenic doped silicon (Si: As) array, 
with a pixel scale of 1.8"/pixel. The slit in the first order 
had a length x width of 57"x3'.'7. We obtained a set of 
16 spectra at each of two default nod positions along the 
slit. The integration time for each individual spectrum 
was 60.95 s, and the total on-source time was 32.5 min. 
In addition, since our target is an extended source, we 
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used the exact same observing mode and obtained spec- 
tra of a nearby sky region for background subtraction. 

2.2.2. Data Reduction 

We started with the BCD images, produced from the 
SSC IRSX Channel-0 Software (version S18.7.0). We first 
combined all the sky BCD images into one master sky 
image and subtracted this master image from each of the 
target BCD images. We then used the IRSCLEAN_MASK 
program, provided by SSC, to construct a super mask 
by combining the rogue pixel mask for the campaign 
IRSX009500 (in which our IRS observations were made) 
and additional bad pixels identified by IRSCLEAN_MASK. 
The super mask was used in the IRSCLEAN program for 
cleaning up bad pixels in the BCD images. The BCD 
images were combined into one final image by using the 
coa2d.pro program. 

We used SPICE, an IRS spectrum extraction software 
provided by SSC, to extract the spectrum of the nebula. 
The width of the extraction box was 15 pixels (27"), 
which covers the 2-D spectrum image well. The spec- 
tra at two nod positions were extracted and average- 
combined into the final source spectrum. 

2.3. WISE Imaging 

Launched on 2009 December 14, the Wide-field In- 
frared Survey Explorer (WISE) mapped the entire sky 
at 3.4, 4.6, 12, and 22 /^m (called Wl, W2, W3, and 
W4 bands, respectively) in 2010 with FWHMs of 6.1", 
6.4", 6.5", and 12. 0" in the four bands, respectively (see 
IWright et alllmOl for details). The WISE all-sky images 
and source catalogue were released in 2012 March. We 
downloaded the WISE image data of the source field from 
the Infrared Processing and Analysis Center (IPAC), and 
found that the cometary nebula was detected in the W3 
and W4 bands but not in the first two bands. The WISE 
observations were made between 2010 Feb. 22-25 and 
the depth of coverage was ~13.7 pixels (corresponding 
to 110 s on-source integration time). 

2.4. Ground-Based Observations 
2.4.1. Optical/ Near-Infrared Imaging 

We also observed the target field in R broad band 
and Hq, narrow band on 2008 June 15, using Inamori- 
Magellan Areal Camera and Spectrograph (IMACS) on 
the 6.5-m Baade Magellan Telescope at Las Campanas 
Observatory in Chile. The long camera (known as f/4) 
of IMACS was used. The detector of the camera consists 
of 8 2048x4096 pixeP SITe ST-002A CCDs. Under the 
f/4 imaging mode, the camera provides a field coverage 
of 15. '4x1 5. '4 and a pixel scale of O'.'lll pixel" ^. The ex- 
posure times of the images were 4 min in R and 18.3 min 
in Hq. The observing conditions were good, with see- 
ing being 0"65 . In a ddition, we observed standard stars 
PG1633+099 (iLand olt 1992^) in R band and LTT7379 
(|Hamuv et al.lll994D in Hq for flux calibration. 

We included near-infrared (NIR) Ks imaging data of 
the target field, obtained on 2006 May 17 as part of our 
IR survey of a few young pulsars (like J1549). The ob- 
servation was also made with the Baade Magellan Tele- 
scope. The NIR camera was Persson's Auxiliary Nas- 
myth Infrared Camera (PANIC : iMartini et al.,.2004il .The 



detector was a Rockwell Hawah 1024x1024 HgCdTe ar- 
ray, having a field of view (FOV) of 2'x2' and a pixel 
scale of 0'.'125 pixel" ^. The total on-source exposure 
time was 22.5 min. During the exposure, the telescope 
was dithered in a 3x3 grid with offsets of 10" to obtain 
a measurement of the sky background. The observing 
conditions were good, with 0"4 seeing in Kg. 

2.4.2. Data Reduction 

We used the IRAF data analysis package for data 
reduction. The images were bias-subtracted and flat- 
fielded. In addition for the NIR data, a sky image was 
made by filtering out stars from each set of dithered im- 
ages in one observation. The sky image was subtracted 
from the set of images, and then the sky-subtracted im- 
ages were shifted and combined into one final image of 
the target field. 

3. ANALYSIS 

3.1. Astrometry 

The pointing of the IRAC frames is typically accurate 
to 0'.'5. In order to determine more accurately the po- 
sition of the cometary nebula and locate J1549 on the 
images, we astrometrically calibrated the IRAC 8.0 /im 
image by matching stars detected in this image to Two 
Micron All-Sky Survey (2MASS; ISkrutskie et "all[200l 
stars. For this calibration, 190 2MASS stars were used. 
The nominal uncertainty of the calibrated image is dom- 
inated by the 2MASS systematic uncertainty (~0"15, 
with respect to the International Celestial Reference Sys- 
tem) . The other images were also positionally calibrated 
by matching them to the 8.0 /im image. 

3.2. Morphology 

We studied the cometary structure by semi- 
quantitatively determining its shape. The 8.0 /im 
image was used for the detailed study because it has the 
best resolution among the Spitzer and WISE detections. 
We first removed all point sources near the cometary 
structure. Using the software package APEX (multi-frame 
extraction) provided by SSC, point sources in the field 
were identified and subtracted from each BCD image. 
The residual frames were then combined into a PBCD 
mosaic, which is shown in Figure [5] 

We then manually selected small regions along the 
three parts of the nebula: the two wings and the head re- 
gion. Narrow box regions set parallel to the three nebula 
parts were defined, with each box having a length x width 
of 5x1 pixeP (6"lxl"2). Photons within each box were 
counted and the uncertainty on the photon counts was es- 
timated as the square-root of the total counts. This way 
a brightness curve that was perpendicular to the narrow 
box regions (or to the nebula parts) was obtained. Using 
a Gaussian function, we performed least-squares fitting 
to the brightness curve and thus determined the location 
of a brightness peak point. The uncertainty for each loca- 
tion was also estimated from the fitting. In total 19 such 
brightness peak points were determined (see Figure ^ . 

3.3. Photometry 

We obtained the flux of the head region of the cometary 
nebula at each band. A circular region with a radius of 
21", which well covers the head region, was used as the 
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Fig. 2. — Spitzer IRAC 8.0 /im PBCD image of th e co metary nebula with in-field stars removed. The small circles mark the brightness 
peak points we obtained for the nebula (see Section |3.2| I. The position of the pulsar J1549 is marked by a black cross sign. The dotted 
curve indicates the 3a region in whi ch the associated object is considered to be located (obtained from fitting the analytic solution for 
bow shocks provided bv I Wilkinll 1999 ) . with the large plus sign marking the best-fit position and dashed curve being the best-fit. T he long 
dashed curve indicate the same region obtained from fitting the analytic solution for pulsar bow shocks provided bv lVigelius et all 1.2007 ). 
with the small plus sign marking the best-fit position and dot-dashed curve being the best-fit. 



photometry aperture. A box region of 80" x 50" near the 
head was used for estimating the sky background bright- 
ness. For photometry at 8.0 /xm, aperture correction was 
required because IRAC imaging uses point sources for 
flux cahbration. We used the function form provided by 
SSC to calculate the correction factor and found a value 
of 0.83 for a radius of 21". The flux value at 8.0 /Ltm 
is the measured flux value times the correction factor. 
We also derived the flux upper limit for the 4.5 fim non- 
detection, because the deep image provides an additional 
constraint on MIR emission from the cometary nebula. 

Following the guidelines given for aperture photometry 
on the WISE atlas imagefl WISE W3 and W4 fluxes 
were also obtained. The broad-band flux measurements 
and 4.5 /im upper limit are given in Table 1. 

4. RESULT 



4.1. Location of the Associated Object 

The shape of the IR cometary nebula indicates that it 
is likely a typical bow shock. The location of the asso- 
ciated object that drives the IR nebula can thus be esti- 
mated by fitting the nebula with a standard bow-shock 
shape, since the distance from an associated object to 
the apex of the associated bow shock (so-called "stand- 
off distance" ) sets the length scale of the bow shock. W e 
considered an analytic solution derived bv lWilkinI ()1996D . 

rg sin 9 — Ta -^^3(1 — 0/ tan0), where rg is the distance 
between the bow shock shell and the associated object, 
is the polar angle from the axis of symmetry (i.e., the as- 
sociated object's moving direction), and is the standoff 
distance. With this solution, we fit the positions of the 
brightness peak points obtained in § 13.21 The free param- 
eters in the fitting were x, y location of the associated 
object, r^, and the rotation angle of the axis of the bow 

® http: / /wisc2. ipac.caltech.edu/docs/release/allsky/expsup/sec2_3f. html 
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Fig. 3. — Dereddened spectrum of the head region of the cometary nebula, which consists of the Spitzer broad-band 4.5 fira fiux upper 
hmit, 8.0, 24, and 70 fim fluxes (diamonds) and IRS 7.5-14.4 /im spectrum (sohd histogram), and WISE 12 and 22 /im fluxes (triangles). 
Dust emission from two temperature components (250 K and 55 K, both dot-da shed curves; t heir s um, long-dashed curve) can generally 
describe the spectrum. The model spectrum of irradiated dust emission given bv lDraine He Lil l|2007l ') is plotted as the red solid curve. 



shock in the plane of the sky. To simplify the fitting, we 
assumed that the axis of symmetry of the bow shock lies 
nearly in the plane of the sky (note that an inclination 
angle such as 60° will result in a smaller projected dis- 
tance between the associated object and the 'apparent' 
apex of the bow shock) . 

The obtained (Str) region in which the associated ob- 
ject is located is indicated as a dotted contour in Fig- 
ure [21 It is a small region very close to the head of 
the cometary nebula, and approximately 15" away from 
J1549 (whose position is R.A. = 15h49"^21!15, Deck = 
-48°48' 37"4, equinox J2000 0, and has an uncertainty 
of iD'Amico et al.l[l998l ). Because the data points 
of the brightness peaks do not provide tight constraints 
along the axis of symmetry, the region is more extended 
along this direction. The best-fit location and bow shock 
solution are shown as a plus sign and a dashed curve, re- 
spectively, in Figure [2] As can be seen, the standard 
bow shock shape generally fits the data points, although 
the tail of the nebula is wider, deviating away from the 
model. From the best fit (x^=14 for 15 degrees of free- 
dom), we found ra — 2'.' 9. 

4.2. Mid-IR Spectrum 

In Figure [31 we show the broad-band and Spitzer IRS 
spectrum of the head region of the cometa ry nebula, 
which is dereddened with E{B — V) — l.l (Sc hlegel et alJ 
1991 using th e interstellar extinction laws given by 
Indebetouw et al. (2005; ) for w avelengths <8 /Ltm and by 
Weinga rtner fc Drainei ()200l D for wavelengths >8 /im. 
The Spitzer /IRAC 4.5 /zm flux upper limit is also in- 



cluded. The IRS spectrum was flux calibrated by equal- 
ing the average value of its 7.5-8.5 ^m region to the 
IRAC 8.0 /im broad-band flux value. A few emission 
lines were detected in the IRS spectrum, which are poly- 
cyclic aromatic hydrocarbon (PAH) 11.3, 8.6, and likely 

7.6 /im. These PAH lines are commonly seen in many 
emission nebulae and arise fr om vibrational modes of 
PAH in dust fe.g.. lDrain^l2003D . Two Ha rotational lines 
were also detected, although weakly. They are H2S(3) 

9.7 iim and H2S(2) 12.3 //m. In general, H2 lines are 
associated with photodisso ciation regions (PDRs; e.g., 
iHollenbach fc T ielcns"1997^ or shock interaction regions 
(e.g.. lShinn et al. .2009) . We note that the broad dip at 
~9.7 /im is not the well-known silicate absorption fea- 
ture. Rather, its appearance is caused by the nearby 
strong PAH emission features praine fc Lill2007[ ). The 
total Galactic reddening along the line of sight to the 
source is only approximately E(B — V) = 1.1 (for 
/=330?5, &^4?3: [SchIe"gel et alJ[T998[) . If the dip is caused 
by interstellar reddening, E{B — V) should be > 4, which 
is too large to b e consistent with tho se estimated from 
the all-sky dust (Schlegel et al."1998') or neutral hydro- 
gen density (Dickcv & Lockman 1990) maps. Also, as 
shown below in this section, the local column density of 
the cometary nebula is much lower than that implied by 
the reddening E{B - V) > 4. 

We first tested usin g a simple dust grain model 
(jPwek fc Werneii I198H ) to reproduce the spectrum we 
obtained. This model was developed to explain the ob- 
served IR emission from supernova remnants (SNRs) , in 
which dust grains are heated both by collisions with am- 
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Fig. 4. — Optical R {left panel) and Ha {right panel) images of the field of the IR cometary nebula. In each panel, the region in which the 
associated source should be located is indicated by the solid curve, while the position of J1549 is indicated by the plus sign. Photometry 
of the field stars (marked by circles in the left panel) was performed, and star a and b (marked by cross signs) was slightly saturated in the 
R and Ks image, respectively. 



bient thermal gas and by radiation flux, and could work 
for the IR nebula in our case. Considering the general 
case that the grains' absorption/emission efficiency de- 
creases as A~^, we found that two components with tem- 
peratures of Tfi = 250 K and T; = 55 K (assuming a sim- 
plified, uniform grain size of radius 0.1 /xm), can generally 
describe the spectrum (Figure [3]). The dust masses thus 
simply estimated are l.OxlO^^d^p^, g and 1.9xl0^^d^p^ g, 
respectively, where dkpc is the source distance in units of 
kpc. 

To better reproduce the IRS spectrum, a detailed 
model including PAH particles is needed. We considered 
the m odel spectra derived and provided bv IDraiiie 
(gOOl), ahhough we note that in their model, dust grains 
are heated by a field of starlight. According to their cal- 
culations, for the case of dust grains being heated by a 
single radiation intensity, the flux ratio between 24 and 
70 /Ltm is mainly sensitive to radiation with energy den- 
sity per unit frequency C/C/isrf, where J7isrf is t he inter- 
stellar radiation field (ISRF) estimated by Mathis^t al.| 
(|1983[ ) for the solar neighborhood and [/ is a dimcn- 
sionless scaling factor. We thus found that J7 = 0.1 
is required for our data points. However, the U = 0.1 
model spectra have much stronger PAH features than 
those in our IRS spectrum (even with the lowest PAH 
mass fraction qpAn = 0.47%; qpAH — 4.6% is found to 
be applicable to the dust in the Milky Way). In order 
to have a good fit to the IRS spectrum, we found that 
the dust should be heated by a single radiation intensity, 
U = 12, plus by a power-law distribution of starlight in- 
tensities ranging from U = 12-10^ (where qpAn — 4.6%, 
the power-law index a is fixed, a — —2, and the frac- 
tion of the dust mass that is exposed to the distribution 
of starlight intens ities is 7 = 0.16; cf. equation (23) in 
IDraine fc Lill2007l) . This model spectrum is displayed in 



Figure [3l which generally describes the PAH features in 
our spectrum. As a result, the dust mass Mdust can be 
estimated from the total va lue of the 24, 70, a nd 160 /im 
fluxes (cf. equation (34) inl Praine fc Lil[2007t ). Since no 
flux measurement at 160 fim was obtained, Mdust should 
be > 1.2 X lO^^d^p^ g and ~ 1.9 x lO^^d^^^ g when as- 
suming the flux at 160 /xm is equal to that at 70/im. The 
value range is consistent with that derived above from 
the simple dust grain model. 

We further estimated the properties of the head region 
of the bow shock. Consideri ng again the dust size o f 0.1 
/zm (and density of 3 g cm~'^: lDwek fc Wernedll98"ll ). the 
total luminosities in the region were 5.7x lO^^d^p^. and 
l.lxlO^^d^p^ erg s~^ for the 250 K and 55 K dust grains, 
respectively. Simply assuming that the head region is a 
21" radius (or O.ldkpc pc) sphere, the number density 
(dominated by the cold dust grains) can be estimated to 
be ^^90 cm^'^, where a mass ratio of 1 to 100 between 
the dust to gas is used (jPraine et al.ll2007l ). 

4.3. Ha Non-Detection 

In order to search for possible line emission from the 
bow shock, we made both R and imaging observa- 
tions of the target field with the Magellan telescopes; 
the wide-band image was used to serve as the continuum 
flux. We did not detect any nebula in our observa- 
tions. Both optical R- and Ha-ha,nd images of the target 
held are shown in Figure HI The i?-band detection limit 
was ~23 mag, and the Ha flux upper limit (3cr) was 
2.9x10"^^ erg s~^ cm~^arcsec~^. 

5. DISCUSSION 

Provided with the results we have obtained from the 
above analyses, we discuss the possible origin of this bow 
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shock by first considering a stellar object and then J1549 
as the associated star. 
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above fe.g.. iCox et al.|[2012[ ). the required ambient den- 
sity is riH — 15; reasonable adjustments to the param- 
eters can easily lead to values more consistent with the 
ISM density. These calculations thus show that it is pos- 
sible to have either a massive star or a giant star with 
mass ejection as the object that drives the bow shock. 

However we evaluated the possible association of the 
bow shock with any in-field stars and obtained negative 
results. First there were no stars detected in the associ- 
ated object region in our optical images down to i? ~ 23 
(see FigureSl). We also analyzed the nearby stars by con- 
structing a color-magnitude diagram of Kg versus R — Kg 
(Figure [S]) . As can be seen, the stars are consistent with 
being low-mass main sequence stars (or white dwarfs) at 
distances of <2.5 kpc. Massive or giant stars would have 
to be at distances of ~6 kpc and ^^10 kpc, respectively, 
and approximately 0.5-0.8 kpc away from the Galactic 
plane, a place not likely to form bow shocks due to the 
low density of the ISM. The detection of H2 lines in the 
IRS spectrum suggests that the associated object should 
probably be < 0.1 kpc away from the Galatic plane, since 
molecular gas generally inhabits within the mid-plane of 
the Galactic disk (Cox 2005). 



Fig. 5. — Ks versus R — Ks diagram of the in-field stars near 
the bow shock. Data points marked by squares are those stars 
highhghted by circles in the left panel of Figure |4] The total 
Galactic extinction of Ay = 3.4 is indicated by the arrow. The 
two stars marked with arrows are those two nearby stars slightly 
saturated in our R or Kg images (FigureSl, with the arrows indi- 
cating their magnitude or color limit. Color-magnitude curves for 
main-sequence stars, red giants, and white dwarfs are plotted as 
the dashed, dotted, and long-dashed lines. 

5.1. A Stellar Bow Shock? 

As shown by our analysis in § 14. 1[ consideriii g an 
isotropic win d/outflow from an object (jWilkinll 19960 . the 
location of the associated object would likely be within 
a small area near the head region of the bow shock, 
which is more than 15" away from PSR J1549. Such 
IR bow shocks can be driven by either m a ssive stars 
with a strong wind (|van Buren et all 119901: iPeri et al.l 
2OI2D or giant stars ejecting mass (jMartin et al.l 120071 : 
Cox et al.ll2012[ ). as long as the associated stars are mov- 
ing supersonically in the ISM. The stand-off distance of 
such a bow shock is given by Va — {mwVw / ^TTpV^Y^"^ 
re.g.. lWilkinl[T99fih . where rh^ is the mass loss rate, Vw 
is the wind velocity, and p is the density of the ISM. 
For a massive 0/B star with typical values of m^, = 



10-^ Mp 



YL. 



1000 km s-\ and K = 30 km s" 



(see, e.g.. iPeri et al.l [2"012l ). the required ambient ISM 
hydrogen density to produce the observed bow shock is 
7iH ^ 1500(mt„_7K.3/K3o^kpc) cm"3, where m„_7 is 



the mass loss rate in units of 10^ Mq yr~ , Vw^ is the 
wind velocity in units of 1000 km s~ , and y*3o is the 
star's space velocity in units of 30 km s^^. Here we have 
used Ta = 2'.'9d = 4 x lO^^dkpc cm and p — pmnnH with 
Ph = 1.4 {ph is the mean nucleus number per hydrogen 
atom and mH is the mass of the hydrogen nucleus) . Us- 
ing the low end value of mu,_7 — 0.1 and dkpc = 10 (see 
below for a discussion of the distance), along with the 
other mass-loss parameters given above, we find tlh ~ 1, 
broadly consistent with values expected for the ISM den- 
sity. Alternatively, for a giant star with mass ejection, 
Vw — ^0 km s~^ and other parameters similar to those 



5.2. PSR JI549 as the Associated Object? 
Given their high space velocities (>100 km s~^^ 



and 



strong winds, pulsars can naturally drive bow shocks in 
the ISM. Detections of bow shocks around pulsars at mul- 
tiple w avelength ranges have provided suIBcient ev idence 
for it (IKasDi et al.l I2OO6I : IGaensler fc Sland l2006h . The 
existence of IR bow shocks, arising due to irradiation or 
coUisional heating, can be expected. The pulsar J1549 
has a spin-down en ergy oi E = 2.3 x 10"^ "^ erg s~^ and a 
distance of 1.5 kpc (iD'Amico et al.l[l998| ). implying that 
if it were associated with the bow shock, approximately 
1.6% of the energy is used to illuminate the dust, which 
is a reasonable fraction to consider (note the total lumi- 
nosity from the head of the bow shock is ^ 3.8 x 10^^ 
erg s~^ at the distance of 1.5 kpc). On the other hand, 
the pulsar's kinematic energy can also contribute to dust 
heating. Considering a radius r — 0.15 pc (the head re- 
gion we have defined when the distance is 1.5 kpc) inter- 
action area, the energy received by the ambient medium 
is il/2pV'^)-V ■{nr'^) « 6.4 x lO^^njjFaoo erg s~\ where 
riH — 1 cm~'^ and V = 200 km s~^, the pulsar's ve- 
locity, are used. Therefore there is sufficient energy to 
illuminate the bow shock. 

The obvious challenge to this interpretation is J1549's 
location, which is at least 15" (or 0.11 pc if at the pulsar's 
distance) away from the determined area for the associ- 
ated object of the bow shock. However, we note that a 
pulsar wind can be highly anisotropic and for peculiar 
cases , long jet-like outflows have been detected fro m pul- 
sars ((Johnson fc WandlMol: [Pe Luca et al.ll20lH) . Cal- 
culations for such anisotropic wind cases to produ ce bow 
shocks have been considered (e.g., lWilkinlT2000D . One 
feature that has been revealed is the off-axis location of 
an associated star if the star's wind points away from the 
direction of the stellar motion, for which the star's loca- 
tion would otherwise be determined to lie in the apparent 
axis of symme t ry when simply fitting the stand-off region 
(|Wilkinll200Cil) . iVigelius et all pOOl did further detailed 
simulations by considering the anisotropy of pulsar winds 
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Wang et al. 



for pulsar bow shocks, and have shown that a variety of 
shape changes with respect to the 'standard' shape of 
bow shocks. In particular, the broadening of the tails 
of the bow shocks in their simulations SA and SB (see 
Figure 13 in their paper), when a pulsar runs through a 
wall of the high-density ISM with finite width, is possibly 
observed in our case (see Figure [2]). We thus did test fit- 
ti ng to our cometary ne bula by using equation All given 
in iVigelius et al.l ()2007[) , an analytic solution for jet-like 
outflows (A, Co, and C2 were simply set to be 0, 0, and 5, 
where A sets the angle between the directi ons of the pul- 
sar's wind and its motion; for details see IVigelius et alJ 
I2007D . The best-fit solution (x^ = 13 for 15 degrees of 
freedom) is shown in Figure [31 which is nearly the same 
as that found in Section ITT] The allowed 3(t region for 
the associated object is however substantially elongated 
toward J1549. While the pulsar is still 5.8" away from the 
region, the fitting results suggest that it may be possible 
for the pulsar to drive the bow shock if it had a highly- 
collimated wind and the wind were tilted away from the 
direction of the pulsar's motion. In order to fully explore 
this possibility, detailed hydrodynamic simulations, such 
as those carried out by Vigclius ct al. (2007), are needed. 

In addition, if J1549 did have a highly collimated pul- 
sar wind, the interaction of its wind with the ambient 
medium could be detectable at X-ray energies. X-ray 
imaging should be carried out to search for the interac- 
tion, which would lie from J1549 to the stand-off region 
of the IR nebula, thus verifying the association of the 
pulsar with the nebula. The bow shock has a direction 
towards the Galactic plane (Figure [T|), implying that the 
associated object is moving in the same direction. High- 
spatial resolution radio imaging of the pulsar, to detect 
its proper motion, should also be carried out. If J1549 
is also moving in this direction, its chance of being as- 
sociated with the bow shock would be substantially in- 
creased. 

6. SUMMARY 

We have discovered an IR bow shock around the pulsar 
J1549, and from our follow-up imaging and spectroscopic 
observations and WISE all-sky survey data, a compos- 
ite spectrum of its head region from 8 to 70 /im, which 
includes a Spitzer IRS 7.5-14.4 /im spectrum, was ob- 
tained. In the IRS spectrum, a few PAH and H2 emission 
features were detected. Considering the overall spectrum 
with these features, the emission likely arises from the 
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dust grains in the interaction region heated by collisions 
and/or radiation flux from the associated object. The 
dust mass in the region was estimated to be 2x10^^ g. 
We analyzed the shape of the bow shock by fitting with 
a standard analytic solution for such objects, and the 
derived region in which the associated object would be 
located is 15" away from J1549. This raises the question 
of which object drives this bow shock. No stars were de- 
tected in the associated object region down to i? ~ 23, 
and moreover our studies of the in-field stars indicate 
that none of them are likely bright 0/B stars or red 
giants located near the Galactic plane, which are able 
to drive such bow shocks due to their strong winds or 
mass ejection. Given the considered properties of pulsar 
winds, including a few peculiar pulsar jet-like outflows, 
we suggest that it is possible for J1549 to be associated 
with the IR bow shock if it had a highly collimated wind. 
Hydrodynamic simulations of the interaction of such a 
wind with the ISM can help explore this possibility. In 
addition. X-ray imaging to detect the interaction of its 
wind with the ambient medium and high-spatial resolu- 
tion radio imaging to detect the proper motion of the 
pulsar will help verify their association. 
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TABLE 1 

Observations and flux measurements of the cometary nebula around J 1549 



Observation date Telescope/Instrument Band/ Wavelength Exposure Flux measurement 









(min) 


(mJy) 


2006 May 17 


Magellan/PANIC 


Ks 


22.5 




2007 Sept 12 


Spitzer/IRAC 


4.5 


26.8 


<0.8 






8.0 


26.8 


10.8±1.6 


2008 June 15 


Magellan/IMACS 


R 


4.0 








Ha 


18. 3=^ 




2008 Aug 30 


Spitzer/MIPS 


24 


11.8 


20.9±0.9 






70 


10.5 


82±32 


2009 Sept 14 


Spitzer/IRS 


7.5-14.4 


32.5 




2010 Feb 22-25 


WISE 


3.4 


1.8 








4.6 


1.8 








12 


1.8 


11.4±0.6 






22 


1.8 


21.1±1.6 



The surface flux upper limit in Ha band is given in § 14.31 
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